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ABSTRACT 

Fusion reaction rates for DT, DD and DHe^ reactions 
have been calculated for non-Maxwe Ilian distribution of 
velocities. In particular, <ctv> values at various kinetic 
temperatures have been calculated for a two - tenqperature 
pseudo Maxwellian distribution, which may approximately 
represent the velocities in a magnetic jnirror reactor. 

It is foimd that for the same kinetic ten^erature, 
a two - temperature distribution gives enhanced reaction 
rates between 5 to 15 keV for DT reaction and between 20 to 
40 keV for DD and DHe^ reactions. These are ti\e temperature 
ranges of maximum practical interest. Calculations at other 
temperatures are also reported. It is found that at very 
high temperatures, the two - temperature distribution gives 
reduced reaction rates in contrast to the favourable beha- 
viour at lower temperatures. 



6 


CONTENTS 


CERTIFICATE 2 
ACKNOWLEDGEMENT ^ 
ABSTRACT 5 
NOMENCLATURE 9 
LIST OF FIGURES 

LIST OF TABLES 13 


CHAPTER 1 
INTRODUCTIOJI 


CHAPTER 2 

ANALYSIS FOR REACTION RATE CALCUIoATIONS 1 

2.1 The Basic Equ^'itions for Binary Reaction rates, 16 

2.2 Assumptions and Their Validity for Fusion Reactions, 19 

2.3 <av> for Maxwellian Distributions, 22 

2.3.1 Three Dimensional Maxwellian Distribution, 22 

2.3.2 One and Two Dimensional Maxwellian Distributions, 25 

2.4 Transformation to CM and Relative Velocity Coordinates, 26 

2.5 A Two— Temperature Pseudo-Maxwe Ilian Distribution, 28 

2.6 <av> for the PseudoHyiaxwe Ilian Distribution, 30 

2.7 Basis of Comparison for Different <av>*s, 34 

2.8 Transformation to Convenient ^Dimensionless) Parameters, 35 


2,9 Summary, 37 



7 


CHAPTER 3 

NUMERICAL COMPUTATIONS AND RESULTS 39 

3.1 Fusion Cross Sections, 39 

3.1.1 Comparison of Experimental and Theoretical Data from 
Refs. [lj, L7J and i,8J, 39 

3.1.2 Interpolation Between Experimental Values, ^6 

3.2 Some Remarks on Computer Calculations, ^6 

3.3 Results and Discussions, 49 

3.3.1 <av> for 1-D, 2-D and 3-L> Maxwellian Distribution, 49 

3.3.2 <av> fora Two-Temperature Pseudo-^'Iaxv;ellian Distribution, 54 

CHAPTER 4 

SUMMARY AND SUGGESTIONS 66 

4.1 Summary, 66 

4.2 Conclusions, 67 

4.3 Suggestions, 67 


APPENDIX A 

A1 Proof of d^v^ ( 1^2 = 

A2 Proof of ’^B^i “ ^*^i 


72 



APPENDIX B 


VALUES OF and a ^ USED FOR CALOULATIONS 

DHe^ 

REPORTED IN THIS V/ORK j^FROM REF [ij 


APPSJMDIX C 

TABLE OF SOME CONSTANTS AND COWERS ION FACTORS USED IN THIS 
WORK 


REFERENCES 



9 


N0I€1HCLATURE 


V 


V 




V 


“A 

mg 

M 


a(v) 

T 


^A ^^1^ 




Relative speed of the two species 

Relative velocity of the tv/o species 

Velocity of the centre of mass 

Nijmoer density of species A per unit volume 

Number density of species B per unit volume 

Mass of the species A 

Mass of the species B 

Total mass 

Reduced mass 

Microscopic cross section for the binary reaction 

Thermod 5 mamic Teaiperature in energy 

Radial temperature 

Axial temperature 

Temperature ratio parameter 

Normalized Maxwellian distribution function of 
speed for species A 

Normalized Maxwellian distribution function of 
speed for species B 



Normalized Maxwellian distribution function 
of velocity for species A 


Normalized Maxwellian distribution function 
of velocity for species B 



11 


LIST OF FIGURES 


FIG. No. 

TITLE 

PAGE__, 




la 

Cylindrical coordinates for d^V-integration 



in Eq, 2.24. 

32 

1b 

Spherical coordinates for d^v— integration 



in Eq. 2.24, 

32 

3.1 

Least squares approximation for 0 ^,^ 

44 

3.2 

Outline of the computer program used to 



evaluate <av>. 

47 

3.3 

Variation of <av> with ^ for DT and DD 

reactions at temperatures T,;=5 and 11^=30 



keV. 

55 

3.4 

Variation of <av> with P for DL and DHe^ 
reactions at temperatures Tj^=20 and T|^=20 



keV. 

56 

3.5 

Variation of <ov> with p for DB and DHe^ 
reactions at temperatures Tj^=-50 and Tj ^=30 



keV. 

57 

3.6 

Variation of <ov> with 3 for DT, DD and 

DHe reactions at temperatures Tj^=10, 



# 

> 

CD 

O 

O 

CM 

11 

O 

o 

CM 

II 

E-i 

58 



12 


Fig, No. PAGE No, 

3.7 Variation of <av> with ^ for DT and DHe^ 

reactions at te^nperatures T|^=5C and Tj^=200 


keV. 59 

3.8 Variation of vs, Tj^ with /D as a 

parameter. 62 

5.9 Variation of vs, Tj^ with ^ as a 

parameter. 63 

3.10 Variation of <av> •:? vs, w'ith i3 as a 

DHe^ ^ 

parameter. 64 

3.11 Variation of vs. with ^ as a 

parameter, 65 



13 


LIST Or' Ti.BLES 

TABLE N O. PAGE N o. 

3.1 Comparison of experimentally measured 

values of from references I'lj and 

LTj V 40 

3.2 Comparison of experimentally measured 

values of from references [Ij 

L7J. 41 

3.3 Comparison of experimentally measured 

values of o , from references [ij 
DKe^ 

and 1.7 j 42 

3.4 <0v>j^^ for 1-D, 2-D and 3-D Maxwellian 

velocity distribution. 50 

3.5 Maxrv;e Ilian 

velocity distribution , 51 

3.6 <av> for 1-D, 2-D and 3-D Maxv/e Ilian 

DHe^ 

velocity distribution. 52 



Chapter 1 


INTRODUCTION 

As is v/ell knovm, in both the lightest and the 
heaviest nuclei, the mean binding energy per nucleon is less 
than that in the intermediate ones. Consequently, nuclear 
energy is released in the fission of heavy elements and in 
the fusion reactions of the lightest ones. The fusion reac- 
tions of practical interest at present in the nuclear fusion 
program are the following : 


D + T — 

— r He^ (3.5 MeV) 

+ n (14.1 MeV) 

(1.1) 

D + D - 

^ T (1 MeV) + 

p (3 MeV) 

(1.2a) 

D + D — 

— > He^ (0.8 MeV) 

+ n (2.5 MeV) 

(1.2b) 

3 

D + He - 

He^ (3.7 MeV) 

+ p (14.6 MeV) 

(1.3) 


For the above reactions to occur, the two reactant 
particles must overcome the Coulomb barrier. In thermo- 
nuclear fusion this is achieved by raising the temperature 
of the reacting species. For useful reaction rates, the 
temperatures required are 5 to 10 keV for the DT reaction, 
and 20 to 40 keV for the DD and DHe^ reactions [l]. These 

* O 

are extremely high temperatures, of the order of 10 K, and 
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most of the difficulties of achieving nuclear fusion are 
related to this. In contrast, we note that a fission reac- 
tion between a neutron and, say, a uranium nucleus, can 
proceed at room temperature, because neutron is not a charged 
particle. 

The reacting mixtures in Eqs. 1.1 to 1.3 become fully 

O 

ionized plasma at temperatures of the order of 10 K. Using 
the cross sections for the reactions under consideration, the 
reaction rates can be calculated if the velocity distributions 
of the reacting species are known. The reaction rates are 
generally calculated by assuming Maxwellian distribution of 
velocities of the two reacting species [lj. 

Our purpose in this work is to calculate the reaction 
rates of the reactions listed in Eqs. 1,1 to 1.3 if the 
velocity distributions deviate from the Maxwellian. This can 
be considered of some practical interest in view of the fact 
that velocity distributions in fusion reactors may be generally 
considered as different from Maxwellian, even though close to 
it. 

For this purpose, v/e analyze in detail the methods of 
reaction rate calculations in Chapter 2, and develop appropriate 
equations for calculating the reaction rates under non-i-Iaxwe Ilian 
velocity distributions. Detailed calculations are carried out 
for a two - temperature pseudo - Maxv;ellian distribution. The 
results are reported and discussed in Chapter 3* 



Chapter 2 


ANALYSIS FOR REACTION RATE CALCULATIONS 

2.1 THE BASIC EQUATIONS FOR BINARY REACTION RATES 

We consider a system containing any mjimber of species 
of particles, and direct our attention on two particular 
species A and B, and infact on a particular binary reaction 
between these two species. Let n^^ and n^ be the number 
densities per unit volume of the two species, respectively, 
at a particular point under consideration, and f ( v^ ) and 
fg ( ) be their normalized velocity distribution functions 

at the same point, i.e,, 

^ ~^2 ^ ~ ( 2 . 1 ) 

Here d v^ and d V 2 are infinitesimal volume elements in the 
velocity spaces. For examjjle, in cartesian coordinates, 
d^v^ = dv^^ '^"'^yl * ^^ 2 l * Choosing v^^ for the velocity space 
of species A, and"!?^ for that of species B is for the conve- 
nience of subsequent discussion, in which one wants to 
distinguish between the velocities of two particles even if 
the species A and B are taken to be the same. 

Let a(v) be the microscopic cross section for the 
binary reaction mder consideration, where 
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V 



( 2 . 2 ) 


If one considers, as a target, a single particle of species 
A moving with velocity v^ , and n^fg (v^) d^V 2 j ■V 2 -v^ | = 
ngfg CV 2 ) d V 2 .V as an incident beam of species B, the result- 
ing removal rate from the beam will be 


V a(v) . 


(2.3) 


7 

Multiplying Eq, 2,3 by n^f^ (v^) d-^v^ , the number of particles 
of species A per mit volume with velocities betv/een ^ and 
v^ + dv^, and integrating over the respective velocity spdces, 
we obtain 


R = 


^A^B 


J 


d^v^ d^V2 


V G (v) f 


A 


(v^) 


"B 


(V2) 


(2.4) 


where R is the number of interactions per unit volume per unit 
time of the desired type, and v is as in Eq, 2.2 

An alternative way of obtaining (or looking at) Eq. 2,4 

^ 7 

is to consider ^^^^4 ("v-i) d-^v^.v as the relative path length 

*»■ ^ 

travelled by particles of species A (between v^ and v^ + dv^) 
in a unit spatial volume per unit time, and to consider 
o(v) Hgfg (^) d^V 2 as the macroscopic cross section due to 
particles of species B betv/een ^ and v^ + d^. Since the 
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latter gives the probability of interaction per imit (relative) 
path length travelled, multiplying the two and integrating over 
the velocity spaces yields Eq, 2,4, 

Equation 2.4 for the volumetric reaction rate can be 
written as 


R = <ov> , 


(2,4a) 


where 


<ov> 




V a(v) 


.fA(v;) ^3^-2) 


(2.4b) 


If A and B refer to the same species, ^^^(vxi) and fg(v 2 ) 
will be replaced by the same distributions, f (^) and f(^) , 
and Eq, 2,4a vi/ill read 


R = ^ n^ <0v> , 


(2.4c) 


where n is the total number of particles per unit volume of 
the single species \inder consideration. The appearance of the 
factor 1/2 in Eq. 2.4c is easily explained on the basis of 
every interaction having been counted twice when considering 
a single species of total number density n [l,6]. Alternatively, 
one can first take n^ = n^ = n/2, and consider the interaction 
between these tv/o halves leading to R^ = (l/4) n^ <av>. One 
can then add to it successive contributions from within each 
half and so on, finally leading to Eq, 2,4c. 
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2.2 ASSUMPTIONS AW) THEIR VALIDITY FOR FUSION REACTIONS 

In deriving Eq, 2,4 for the volumetric reaction rate, 

R, various assumptions have been implicitly made. Some of 
the (more obvious) assumptions, and their validity for fusion 
reaction rate calculations, are briefly discussed belov;, 

(1) Eq, 2,4 assumes that the particles of species A 
and B can be considered as classical particles, i,e, as loca- 
lized wave packets v/hose extensions are small compared to the 
average interparticle distance. This will be realized if the 
average de Broglie wavelength, is much less than the average 
interparticle distance, n [12J s 

4d » ( 2 . 5 ) 

^ “mv » (2,58) 

where n is the particle density, h Planck’s constant, m the 

mass of the particles (ions in a fusion reaction) , and v is an 

average speed, say, the rms speed of the particles. If the 

condition 2.5 is not satisfied, the linear dependence of R on 

n^ and n^ (separately) in Eq. 2,4 will breakdown, even if 

suitable definitions of other variables can be arrived at. For 

a fusion reaction at say 10 keV, it is easily seen that 

10 m which is sufficiently less than the interparticle 

31 —3 

distance even at n = 10 n a typical particle density for 
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inertial fusion. (Some of the constants needed for these 
estimations are given in Appendix C ) . 

(2) The range of the forces leading to the reaction 
under consideration in Eq. 2.4 should he sufficiently less 
than the interparticle distances. The best way to check it 
is % 


,1/2 ^1/3 


n 


<< 


(2.6) 


^00 p 

For fusion reactions, the highest value of a is 5x10 m , 
and the condition 2,6 is seen to be satisfied even in the 
worst case of inertial fusion. It may be mentioned that if 
condition 2,6 is not satisfied, again the linear dependence 
of R on n^ and n^ v;ill breakdown, i.e, some other term(s) in 
the rhs of Eq, 2,4 will also have some dependence on n^ and 

ng. 


(3) For Eq, 2,4 to be valid, the velocity of a 
particle should not be correlated to its position with respect 
to another particle. For this to be realized, there should 
be no long-range interparticle forces. This obviously is not 
valid for plasmas. However, if the average potential energy, 
<PE>, in the system per unit volume is sufficiently less than 
the average kinetic energy, <iQS>, per unit volume, the effect 
of correlations on the reaction rate car be ignored [6,4], 

This is equivalent to the condition that the Debye length, kg, 
be large compared to the interparticle distance [2,3J : 
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» 1, (2.7) 

where 

e„ T 1/2 

^ ^ » ^2. 7a) 

n e 

= permittivity of free space, 

T = temperatur*e, in J,and 

e = electronic charge. 

It is seen that for inertial fusion (nw-10^"* T 10 keV) , 

1 /"^ 

Xjj n ' 1, Hence the condition 2.7 is not satisfied, and 

correlation corrections may be required. However, all the 
calculations in this work re3.ai;e to magnetic fusion, in parti- 
cular, to velocity distributions appropriate to mirror systems. 
For magnetic fusion, the condition 2.7 can be easily seen to 
be satisfied, 

(4) It should be mentioned that Eq, 2,4 gives only 

an expected reaction rates. For magnetic fusion smallest 

volume aV and time 5t of interest are of the order of 10~^ 

and 10 s, respectively. For inertial fusion, the same may 
—IP! —11 

be ' — r 10 m^ and v^'% 10 s. At T = 10 keV and densities 

21 — '51 —'5 

of the order of 10 m and 10 m respectively, it can be 

Q 

seen that RAV5t » — ^ 10 , Hence the expected rate ’-jill be close 
to the actual rate, fluctuations in RA'Vdt being of the order of 
sliWVdi [6j. 
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(s) Relativistic corrections are completely negli- 
gible, average ion velocities in fusion plasmas being of the 
6 —1 

order of 10 m s . The effect of external forces on fusion 
cross sections is ignored (Hov/ever, see the next rem*ark) , 

(6) Polarization effects are ignored, i.e,, in the 
meas'urement or calculation of o(v) , all allowable relative 
spin orientations of the interacting species are considered 
to be equally probable. ( Oppositely polarized ions can enhance 
fusion cross sections by as much as a factor of tvro [11J, The 
external forces, for example, magnetic fields, do effect the 
state of polarization, generally acting as a depolarizing 
mechanism, and thus destroying the reaction enhancing effect 
of spin polarization.) If the state of polarization is clearly 
known, an appropriately weighted ct(v) can be defined and Eq, 
2,4 for reaction rate calculations can still be used, 

2,3 <av> FOR M»/ELLIAE DISTRIBUTIONS 


2.3.1 Three Dimensional Maxwre Ilian Dis'bribution 


The three dimensional normalized Maxwellian distribution 
of velocities for, say, a species A, whose (peculiar) velocities 
are denoted by, say, v^ , is given by; 




■°Ah 

m. 3/2 2T 
— ) e 
27rT 


( 2 . 8 ) 
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from which it follows that the tliree dimensional Maxwellian 
distribution of speeds is 




m. 


3/2 


2 T 


2tc 


rn / 


(2.8a) 


As usual, we use the same S 3 nnbol for the velocity and speed 
distributions, even though the two functions are different. 

In Eqs. 2,6 and 2.8a, ra^ is the mass of the particles (in Kg.) 
and T is the thermodynamic temperature (in J), Normalization 
is easily verified by integrating Eq, 2,8a w.r.t, v^ from 

9 

v^ = o to « , Alternatively, in cartesian coordinates, v^ 

2 2 2 

in Eq. 2,8 can be replaced by v + v + v and integration 

x1 y1 z1 

performed w.r.t. v^^ , v^^ and v^^ from - » to f « . In both 
cases, the result obtained is unity. It should also be pointed 
out that the distribution in Eq, 2,8a is for all solid angles. 

If it is written on the basis of per steradian, the factor 471 
will not appear in the distribution function, but will reappear 
during integration over the entire solid angle. 

In order to calculate <av> (see Eq, 2.4b) for two 
species A and B, both of v/hich are at the same thermodynamic 
temperature T, the relative speed distribution between the two 
species can be i-jritten down using the basic laws of statistical 
mechanics [6j, i,e, the same laws which yield Maxwellian distri- 
bution of individual species to begin with. This is obtained 
by replacing m in Eq, 2.8a by the reduced mass, P- i 
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m^+ifTg 


(2.9) 


^AB 


, U 3/2 

4:1 ( ) e 

2TrT 


2T 


v2 


( 2 . 10 ) 


where v = jv^ “ | is the relative speed as in Eq. 2,2, 

The meaning of i'j:^ 3 ('v) is that dv gives the number 

of particles of species B per unit volume whose relative 
speeds w.r.t, any given particle of species A are between v 
and V + dv. Alternatively, (v) dv can be interpreted 

reciprocally in a similar way. With this interpretation, it 
is then straightforward to rewrite Eq. 2.4b for this particular 
case of two species being at the same thermodynamic tempera- 
ture T [l] s 


<0V> 


3D 


4Tr(^) 


3/2 


dv v^ o(v) e 


1£_ 

2T 


v2 


( 2 . 11 ) 


where 3D denotes 3“dimensional Maxwellian distribution, 

¥e may mention here that the argument leading to Eq. (2,11) 
may appear somewhat subtle. In Section 2.4, we will obtain 
expressions for" <av> for more general (non— equilibrium) velocity 
distributions, Eq. 2.11 will then appear as a special case, ¥e 
may also mention here that Eq, 2,11 is the standard equation 
for calculating <av> for thermonuclear reactions [l,5]. 
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2.3.2 One and Two Dimensional Kaxw elli an Distributions 

Although one and tv/o dimensional Maxwellian distri- 
butions are e:itrenie cases and perhaps not of practical 
interest, ■'..e have calculated <0v> for these distributions 
particularly to see that tiiese limiting results are approached 
smoothly in computer calculations, as well as in analytical 
formulations. It is interesting, however, to note that if 
these limiting distributions could be approached in practice, 
it can lead to sigriif leant reaction enhancing effect in 
cases of practical importance (see Chapter 3 for details). 

As has been remarked above, these limiting cases v;ill 
appear as special cases of the more general expression for 
<av> to be developed in the next section. Our purpose here 
is to indicate that the argument used in Section 2,3.1 for 
three dimensional Maxwellian distribution can also be imme- 
diately extended to obtain <av> for one or two dimensional 
cases. For convenience, we denote the temperature T by 
for 1-dimensional case, and by for 2-dimensional case. 

By denoting the speeds by v,j^ and respectively for the 

1-D and 2-D cases, we have the following speed distributions 
for species A i 

\ 

r . ..2 

^A ^1z 

m. 1/2 “2T 

^ ’ ^ 2 . 12 ) 

[1 -dimensional Max'.vellian distribution 

of speeds] 
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and 


DIa 

^27tT ^ 


2 T^ 


V 


1r 


(2.13) 


[2-dimensional Maxwellian distribution ] 

of speeds 


By using the same argijments as in Section 2.3.1, one can write 
respective ^AB (v) by replacing m^ in Eqs, 2.12 and 2,13 by 
the reduced mass, M, and finally the follo\/ing expressions for 
<av> can be obtained i 


and 


cx 


, u 

<av>.j, = 2 (-=-—) 

27tT 

z o 


11'^ 

dv V a(v) e . (2,14) 


'0 


/ 


"^*^^^20 “ ( ip ) -^ 

r o 


dv v^ C7(v) 


2T^ 


(2.15) 


2.4 TRANSFORMATION TO CM AND RELATIVE 'i/ELOCITY COORDINATES 


In order to calculate <av> for two arbitrary velocity 

distributions f^( v^ ) and fg( 2.4b has to be inte- 

— ^ ^ 

grated over six variables - 3 for v^ and v^ each. As it 
stands, none of these integrations can be performed without 


J ' ■ ■' X 

making explicit use of a(v) , because v = | v^ - v^ 


contains 


all the six variables of integration. However, by changing the 

— ^ ^ ) 

variables of integration from v^ and V2 to V and v (defined 
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below) , five of the integrations can be performed as o(v) 
does not depend on them. Integration over v can then be 
performed analytically or computationally, making explicit 
use of a(v) , 

¥e define ; 


IT 


"'a "^1 




m-B v^ 


and 


(2.16) 


V = 


^2 ““ 


(2.17) 


It is immediately seen that V is the center of mass velocity, 
and V is the relative velocity, whose nagnituae v v/e have used 

—Tt — ^ 

earlier. By changing the variables of integration to V and v, 
Eq, 2.4b can be written as 


<cv> 



^(v^, v„) 

1 d^v d^V V a(v) f^(V,v) fg(V, v) 



where 


5 (v^, 

i(^ t) 


(2.18) 


is the magnitude of the Jacobian 


h (v., 

— 3r — I't should be mentioned here that after the variables 

(? , t ) 


are changed, in writing Eq_, 2,18, we have continued to denote 
the distribution functions by f^ and fg, even though they are 
different functions. Further, the region of integration in 2,18 
should corresponding to the region of integration in the original 
Eq. 2.4a. 
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The Jacobian mentioned above is evaluated in Appendix 
A, and its magnitude is seen to be unity. As a result, 


<av> = 


d^v dpY f-r, (“v, ”?) V a(v) 


(2.19) 


It is seen that in Eq, 2.19, integration w.r.t. d^V can be 
performed in any convenient coordinates. However, for inte- 

3 

gration v.r.t. d v,one of the coordinates must be the magni- 
tude of i.e. V, so that tv^'o of three integrations involved 
can be obtained v^ithout using a(v) . Spherical coordinates 

3 

for d V will, therefore, generally be suitable. 


2.5 A TWO-TEHPERATURE PSEUDO-M]r.VE,LLIAN DISTRIBUTION 

In the next section, we will develop an appropriate 
expression for <av> for a chosen nonH>Iaxwe Ilian distribution 
f-unction of speeds. For this purpose, we choose a two- tempe- 
rature pseudo-ilaxwellian distribution function, which in some 
sense, can be used to characterize the velocity distributions 
in a mirror fusion reactor. This is based on having one 
temperature, T^, characteristic of the velocity distributions 
along the magnetic field lines, and another temperature, T^» 
characterizing the velocities perpendicular to the magnetic 
field lines, A norma j.ized velocity distribution function of 
this kind for, say, a species A whose velocities are denoted 
by , can be written as follows % 
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2T^ 2T^ 

e " , (2.20) 

where 'v >|2 = component of parallel to the field lines, and 

Vir = component of perpendicular to the field lines. 

In Eq, 2,20 if we put, 

= T, and 

2 2 2 

= ■'^1r ''^Iz » 




2 

m. 2 m. 

(_A.) (_A. 

27tT 2t[T 


¥e immediately get the 3-D Maxwellian distribution of Eq. 2.8. 
If we multiply Eq. 2,20 by 27Tv^^dv^^ and integrate from 
v^r = o to oo , we get l-D Maxwellian distribution of veloci- 
ties : 




2nT^ 


1 

2 


) e 


“A^lz 

2 


( 2 . 21 ) 


It is to be noted that the velocity distribution of Eq, 2,21 
above rightly differs from the corresponding speed distribution 
of Eq, 2.8 by a factor of 2. This is because the velocity v:^^ 
will range from -» to + «* , while the speed v^^ VJ'iH range 
from o to at # 

Similarly, if Eq, 2.20 is integrated over v^^ (from-W\, 
to + ) , 2-D Maxwellian distribution of velocities is obtained; 



30 





( ) e 

2711 - 


m 




2 



( 2 , 22 ) 


which again differs from the corresponding speed distribution 
of Eq. 2,13 by a factor of 27iv^^. 

2.6 <av> B'OR THE PSEUDO-I'lAXl.'ELLIAK DISTRIBUTION 

In order to calculate <av> for tv;o species A and B, 
with velocity distribution as in Eq, 2,20, we first write a 
similar expression for f^( ), taking and for species 

B to be the same as that for species A, (If this is not so, 
the method developed in Section 2,4, which is used below, will 
still be applicable. Only the details will be different and 
more involved.) 

In general, in order to change to CM and. relative 
velocity coordinates, Eqs, 2,16 and 2.17 will be solved for 
six components of v^ and v^ in terms of the six components 
of V and v, and the resulting expressions will be utilized 
in fj^(^) and f 3 ('^) to obtain f^Cv,”^} and fgC^'^). Then, 
as pointed out earlier, five of the six integrations in the 
general expression for <av> in Eq, 2,19 can be carried out 
without explicit use of a(v) , 

The same procedure is, of course, applicable for the 

distribution being considered. However, the details are 

simplified for this particular case by noting that Eqs. 2.16 
and 2,17 imply for z-coinponent (see Appendix A) i 
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"'a Az '■ ^2.23) 

where ii is "the reduced mass as defined earlier and M = m.+m-r,. 

A B 

Relations similar to Ec, 2,23 ’v.'ill be valid for other compo- 
nents, As a result we also have, 


+ mg V2^ 


d- V, 


2 


M V 


2 


(2.23a) 


9 2 

Since only v^^ and v^^ occur in the pseudoHl'iaxme Ilian distri- 
bution of Eq, 2.20 under consideration here, it is easy to see, 
with the help of Eqs. 2.23 and 2.23a, that the Eq. 2.19 for 
<ov> reduces to the following for the present case ; 


mg) 


3/2 


<av> 


rz 


^ 2 
Stt-T^T/ 



d^Vd^v V a(v) e 


''V 

2 T. 




2 


2T„ 


2T„ 


(2.24) 


where we have put a subscript rz on <ov> to indicate that this 
applies to the two-temperature distribution of Eq, 2,20. 


As promised earlier, five of the six integrations 
involved in Eq, 2,24 can now be performed. This is done by 

"Z 

writing d-'^V in circular cylindrical coordinates (Fig, la) i 


d\ = dVj, d& dV^ , 


and integrating with the following limits ; 
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FIG, la ; Cylindrical coordinates for d''V-integration in Eq,2.24 
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FIG,1b ; Spherical coordinates for d^v-integration in Eq, 2.24 
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s o to « 

; - 08 to 30 , and 

9- t o to 271 (see Fig, 1a) 

Two of the integrations related to d^v can be performed by 

*5 

writing d v in spherical coordinates (FIG, 1b) ; 

d\ = sin9d9d0- , 

and replacing and in Eq, 2,24 by : 

= V cos 9 , 

= V sin9 , 

with limits i 


9 ; 0 to TT 

0 s o to 271 (see FIG, 1b) 

The resulting expression for <av> ^ is i 

Jl ^ 


<av>. 


rz 


■jr 


3/2 




T xl/2 

z r ^ 
o 


C!(v) 


2 T^ 


e , 


I(v) dv 


(2.25) 


where , 


I(v) = 


iLv2(T^-T^) 

o rp 

z r 


, c3a 


( 2 . 26 ) 
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[a here is related to the original integration w.r.t. 9 by s 
a = cos 9 5 this integration can not be done analytically, but 
it does not require the knowledge of a(v)]. 

If we use (= T) in Eq, 2*25, we immediately 

obtain as in Eq, 2,11, It is also possible to take 

in Eq, 2,25 as T^. . . ^ o, or as o, and obtain 

<av>^D respectively. However, 

these limits is somewhat involved, and the desired 
results are more directly obtained by starting with Eqs, 2,21 
and 2,22, instead of Eq. 2.20, and repeating the steps out- 
'lined in this section for these special cases. 


2.V BASIS OF C0I4PARIS0N FOR DIFFERENT <ov>*s 
\ 

In order xo compare <dv>'s for arbitrarily chosen 
velocity distributions ^ ) and fg ( ), we define 
effective (kinetic) temperatures, and Tg, of species A 
and B, respectively, as follows : 


2 ^A = 5 ®A ^ ’^1 


and 


I Tb - 2 ”^B ^2^^ 


For the two-temperature pseudoHyiaxwe Ilian distributions [of 
the form given by Eq, 2,20], with the same T and T for the 

JL ^ 

two species, these effective temperatures are easily seen to 


55 




I,' 


if' 

h 



be the same. Denoting tMs common effective (kinetic) 
temperatinre by we find by straightforward integration; 


1 

2 



+ 4 T. 


(2.27) 


Since is a measure of the average kinetic energy per 
particle, it can be used as the basis for comparing <ov> 
values for different distributions. Thus while varying 
and T„ for the two-ternperature distribution, we compare <ov> 
values for a fixed at a time. 


2,8 TRANSFORMATION TO CONVENIENT (DIMENSIONLESS) PARAMETERS 

To facilitate integration of Eqs, 2,25 and 2.26 in a 
more elegant form, and to facilitate some of the further 
discussions, we define 



(2.28) 


v;here jj, is the reduced mass. We then define the temperature 
ratio parameter, by ; 



and a dimensionless relative speed, by : 



(2.29) 


( 2 . 30 ) 
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These definitions, Eqs,,2,27 to 2.30, can be introduced into 
Ea,s, 2,25 and 2,26. We first introduce Eqs, 2,2? and 2.29 to 
obtain t 


<av> 


rz 



v^o(v) ICv,p) 


-( 23 + 1 )^ 


dv 


where 


I(v,p) 



(H-2P)(1-P) q. ^2 

— *2 

6 3 T,. 

e .da (2,32) 


The above equations bring out the dependence of on Tj^ 

and 3. If we now introduce Eqs, 2,28 and 2,30 in the last two 
equations, the later can be witten as ; 


<av>, 


rz 


OO 

= f oM I(^%3) e 


_ ^^( 23 + 1 ) 
23 


. dw , 
(2,33) 


where 


I(w,3) 


1 

■/ 


(1+23) (1-3) 
23 


a' 


da 


(2,34) 


It is easily seen that the case 3 = 1 corresponds to <av>^ 
of Eq, 2,11, It is also possible to show, with some involvement, 
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"tha'fc "the limits of Eq, 2«31 as I3 — and p co exist and 
correspond, respectively,, to and of Eqs. 2.14 

and 2,15, To avoid repetition, appropriate expressions for 
these special cases are direcxly given in the next section, 
where we sumixiarize the results of this chapter v/hich are 
actually used in later calculations. 

2,9 SWiARY 


For a two-temperature (pseudoHyiaxwellian) distribu- 
tion of velocities for species A, 


m. yr m, 

f.(^) = (— )^ . 


2nT 


z 


2nT^ 


2 T„ 


2 T„ 


(2.35) 


and a similar distribution, ig( ) , for species B with same 
T and T , v;e have 


<av>. 


/2 (2j3+l)^'^^ 


rz 


71 * 



_ ■^^( 26 + 1 ) 

2B 

3 .e .dw , 

( 2 . 36 ) 


where 



I(w,^) 


0 


(2.37) 
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3 = ^ 2 T 

^ -^z » 2 


T + ~ T 
r ^ 2 » 


3 T, 


= 


K 


JLL 


CO = 


V 


V = 


Vj - 


, and 


iCL = 


m^ + nig 


For 3 = 1, 


<ov>jD = 3 


I Vj- I o(tJ) 


■2 0)^ 


dw , 


for 3=0 


<av> 


1D 


■Jij- 

' n 


a(w) 


(O' 


d . 


^ » 


and for 3 = «> , 


<av>, 


2D 


oe 

/ 


-CO' 


2 Vj^ I “ cr(co) 


(2.38) 


(2.39) 


dco 


(2.40) 



NUMERICAL COMPUTATIONS iiND RESULTS 


3.1 FUSION CROSS SECTIONS 


3.1.1 Comparison of Experi f ner it al an d Theoretical Data from 




Tables 3,1 to 3.3 sho\'/ experimentally measured values 
of and 0 at ce-lected enei-gy points. Here the 

incident particle is deuteron, and the target particle (triton, 
deuteron and He^ respectively) is considered to be at rest in 
the laboratory system. The deuteron energy shovm is, accordingly, 
the energy of the incident deuteron in the laboratory system. 

It is related to the relative velocity, v, of Section 2,9 by 


^ v^ = (3.1) 

where is the deuteron mass. 

It is seen that there is considerable difference in the 
cross sections from the two soinrces. For all calculations in 
this work we have used the cross sections from Ref. [l]. The 
values in Ref. [lj appear to have been arrived at by consulting 
different sources of experimental data, and were therefore 
considered to be better suited for our pijirpose. The detailed 
values are given in iippendix B. It should be noted that the 
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Ilore detailed values are listed in Append .i 3 



Comparison of experimentally measured values of a from references 

DHe^ 








Note ; The values used for all calculations in this work are those of lIJ. 
More detailed values are listed in Appendix B, 
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cross section is ■the ioial cross section for both the 
reactions shown in Eqs. 1.2a and 1.2b. 


On theoretical grounds, the cross section for a nuclear 
fusion reaction can be ’//ritten as [ijs 



(3.2) 


where W = energy of the ineracting ‘particles in center of mass 
system, and 

A and B are constants. 

Eq, 3.2 is particularly valid if W is \;ell below the Coulomb 
barrier between the interacting particZ-es, Theoretical expre- 
ssions of A and B are available ([lj, [sj). Ho’wever, consi- 
dering A and B as consta'nts to be empirically determined, a 
least square fit between Eq. 3.2 and experimental data from 
Ref. [l] for DT reaction is shown in Fig, 3.1. For this piirpose, 
only seven representative experimental values shown in Fig, 3,1 

were used. The imweighted square error, ^ , in the cross sec- 

s 

tion values, given by, 

Eg = S (a - (3.3) 

was minimized, where o is as calculated from Eq. 3.2 at the 
representative energy points, and is the- corresponding 
experimental value. Here it shovild ba mentioned that if A and 


xIO®). 
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B are considered as empirical constants, ¥ in Eq, 3.2 can 
be replaced by deutercn energy, in the laboratory system. 

It is seen in Fig, 3,1 that there is considerable 
difference between theoretical and experimental values, even 
if the constants A and 3 are determined empirically by least 
squares approximation. In this connection, we may also remark; 

1, In Ref. [sj , a theoretical expression for o(ol>) is 
given as follows s 

- /b/w 

o(v;) = a e V . (3.4) 

¥e considered this expression also, but the least squares 
error was more in comparison v/ith that obtained from Eq, 3,2. 

2, It appears that the values given in Ref, [l] 
actually correspond to a least square approximation between 
experimental data and a plot of In (E^^ a(E,^)) vs Ej 3 
according to Eq. 3.2, This plot is linear as can be seen from 
Eq. 3.2, by replacing W by E^^, ¥e tried to check it. The 
values do not fall on a straightline. (There is some confusion 
here, probably due to a misprint, as in Ref, [l], it is stated 
that the plot In (Ej^ c7(Ej^)) vs E^ is a straightline in accor- 
dance with Eq, 3. 2 1 this is obviously not true.) 

In any case, the cross section values given in Ref,[l] 
are most extensive, and are in reasci able agreement with limited 
experimental data consulted during this work. In Ref . [l], <av> 
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values are also given for three dimensioral Maxwellian 
distribution of velocities. Therefore, for the purposes 
of comparison with some of the calculations reported in this 
work, it was thought best to use the cross section data from 
Ref. [lj, as summarized in Appendix B. 

3.1.2 Interpolation betwee n Experimentol 

For integrating Eq, 2.36, aii energy interval of 1 keV 
was used. Since the values given in Appendix B are at inter- 
vals ranging from 1 to 100 keV, an interpolation had to be 
used. We considered both linear interpolation and an inter- 
polation based on Eq, 3.2, with A and B determined using the 
two end values under consideration. The final results reported 
in this woi’k are all based on this later interpolation scheme, 

3.2 SOME REMARKS ON OOlviPUTCR CALCUIJiTIONS 


An outline of the computer program used for calculating 
<av> in accordance with Eqs. 2.36 to 2,40 is shown in Fig, 

X Cmi 

3,2. All integrations were done using Simpson’s rule [l3Js 


x+h 



f (x) dx 


X 


^ lf{x+h) + f(x)j 


(3.4) 


■ 3 : 

with an error term of the order of ir' , This v/as found to be 
adequate when used \ath an energy interval of 1 keV. However, 
for integrating Eq, 2.37, especially for low values of very 
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small intervals of c: were needed to obtain desired accuracy. 

For most calculations 20 to 50 uniform intervals of a over 
the range o to 1 were adequate. But for small values of 
as many as 1000 intervals of a were needed, Ptirther, it was 
found to be the best to merge the exponential in Eq, 2,36 with 
the ex'ponential in Eq, 2,37 for the purposes of integration, 
and to use non-uniform intervals of a. Relatively large 
intervals in the ran^e o to 0.8 and finer intervals in the 
range 0,8 to 1 gave good and fast results. 

The correctness of the calculations \7as verified by 
the follo\.'ing considerations t 

1, The values of given by Eqs. 2.36 and 2,37 for 

j 3 — and ■ ■» eo were calculated and compared with the limit- 
ing results of Eqs, 2,39 dnd 2.40 respectively. Calculations 
were improved until accurate agreement was achieved upto 
3 = 10""^ to jS = 10^, Since these v;ere the worst cases, the 
accuracy of the calculations for other values of p could 
obviously be assumed, 

2. For P = 1, should correspond to 

<av> 3 jj of Eq, 2,38. This was checlced. These values v/ere 
further compared with the <ov> values given in Ref. [l] for 
3-dimensional Maxwellian distribution of velocities. Expected 
agreement was found, as will be reported later. 
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3. If we put a(w) = 1, then Eqs. 2.36 to 2.40 can be 
integrated analytically or through the use of tabulated error 
function for 3 = 1 , o and •«> . The results thus obtained 
agreed with the computed results for these cases. Here we 
may point out that actual limits of integration v/ith respect 
to w in Eqs. 2.36 and 2.33 to 2.40 are not o to » , but 
correspond to deuteron energy values from 7 to 1000 keV, as 
the data available was limitea to xhis range (see Appendix B) , 

3.3 RESULTS AIID DISCUSSIOr'S 

<av> - calculaxions \/er3 doiia for the DT, DB and 
■5 

DHe fusion reactions for the I'!a 2 :v/e Ilian and non-iiaxwe Ilian 
velocity distributions, ^/ith a viev/ to see the effect on <av> 
of a departure of velocity distribution from !laxv;ellian. The 
departures considered were limited to those vhlch could be 
described by a two-temperature pceudo-i''Iaxv;e Ilian distribution 
as discussed in Section 2,5. The values of <av> for different 
distributions are then compared for the same value of Tj^, the 
kinetic temperature, as discussed in Section 2,7. We summarize 
the results in the following subsections. 

3,3.1 <av> FOR 1-D, 2-D and 3-D KAXV J ELLIA.:J Distribution 

The values of <av> available in literature are for 
3-dimensional Maxwellian distributions. These are listed in 
the last columns of Tables 3,4 to 3.6, taken from Ref, [l]. 




Table 3,4 




<av>-p.T^ for 1-D 2-D and 3“D Maxwellian velocity distribution 


I 



32074 



<av> :z for 1— D, 2-D and :)-D Maxwellian velocity distribution 
DHe^ 
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In "bhe columns ;ius't preceding "the Iasi; columns, are listed 
the results of our calculations [Eq. 2,33j. It is seen that 
there is general agreement, Hov;ever, in view of the fact 
that the input data (Appendix B) used for all our calculations 
is also from Ref, [l], better agreement could be expected, 

¥e do not v/ish ~o stress this comparison too much, except to 
note that if the^input data is as in Appendix B, then <av>^j^ 
values should be as calculated here. 

Our main interest here is to see the effect of depar- 
ture from Ma:orellian distribution on <av>. For this piurpose 
we list in the same Tables 3.4 to 3.6, the calculated values 
of and <av>2Q» as calculated in accordance v/ith Eqs, 

2.39 and 2.40. As already noted in Chapter 2, these are 
extreme forms of the t^vo- tempera ture,pseudo-4'Iaxv;e Ilian distri- 
bution as jS — > o or w respectively. This is hardly a situa- 
tion which would ever exist in practice, but establishes the 
upper boimds of the effect on <av> of possible departures 
from Maxwellian distribution. 

As we see from Tables 3.4 to 3.6, <av>^j^ and <o-v>^ 
can be significantly larger than <av>^jj upto a certain tempera- 
ture Tj^, This temperature is about 10 keV for DT reaction and 
about 30 keV for' DD and DHe-'^ reactions. In view of the fact 
that these are the temperatures of practical interest in thermo- 
nuclear fusion, we arrive at the conclusion that 1-D and 2-T) 
Maxwellian distribution of velocities give enhanced reaction 
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rates £or all the taree fusicr reactions. However j at lighter 
temperatures, the reverse eiiect is seen. 

These effects can be e^’pleined. as fjllows. For the same 
kinetic temperature the average rdotive speed is highest 
in 1-D and lowest in 3-1^ Ma-'-v.ellian ois-:ributions. At lower 
values of this contributes positively to <jv> as o increase 
with V. But beyond a certain v.o decreases with v. This is 
apparent fi-om Append!:; B orJy for DT reaction. However, the 
same thing happens for DD and DHe^ reactions if we push the 
data beyond 1000 keV in Appendix B. In any case, as far as 
numerical calculations in this \.-crk are conceinied, a is taken 
to be zero beyond 1000 keV. 

Thus at higher T^, the gain in v actually reduces <av>. 

3.3.2 <ov> For a TV'/o -Temperature Pse udoH !via:gtfe Ilian 
Dist r ibut ion 

Figs. 3.3 to 3.7 plot of E.-i, 2,36 as a function 

of the parameter 3 = T /T . As can be seen from the figures, 

X 

O A 

while 3 is varied, the kinetic temperature, % = (2 '^z'^'^r^ 

is kept fixed. The meaning of this is, that as the radial and 
axial temperatures are varied, the average kinetic energy per 
particle is kept fixed, Tj^ being a measure for the same. This 
has already been discussed in Section 2,7. 

Figs. 3,3 plots the results for the DT reaction at 5 keV 
and for the DD reaction at 30 keV, the temperatures of maximum 




temperatures T|< = 5 and T|<=30keV 


041 




f c 


ni x<'rTA> 


FIG. 3*5 Variation of <crx^> with p for DD and DHe^ reactions 
at temperatures Tj^ = 50 and Tj<=30keV. 



fN| 



FIG. 3-6 Variation of <Oi^> with P for DT, DD and DHe^ reactions 
at temperatures T|^ =10, T|^ = 200 and T|^ = 100keV. 
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interest; in nuclear fusion for these reactions, respectively. 
In the vicinity of the 3~D Maxv/ellian distrihution, enhance- 
ment in the <-!v> for DT reaction at 5 keV is seen to be more 
than 20% (^ 1 *%^ 0,5 to 2,0) , The enhancement for the DD reac- 
tion at 30 keV is, hov;ever, only about ']%, It is seen xhat 
the reaotion rate is enhanced if no matter whioh is 

more. This oan also be confirmed by differentiating 
in Eq, 2,36 partially p. It shovvs an extremum at 3=1, 

as is also apparent fusm the computer resi;lts shown in Figs, 
3.3 to 3.7. 

Figs, 3Ji- and 3.5 show similar data for DD and DHe^ reac 
tions, at some chosen values of Again an enhancement in 

<av> value due to deviation from kaxwe Ilian distribution is 
seen. 

Fig. 3,6 gives the results for DT reaction at 10 keV, 
another temperature of practical interest for nuclear fusion. 
An enhancement of the order of 105^ in the vicinity of 3=1 is 

3 

noted. In this figure, we also give the results for DKe 
reaction at 100 keV and DD reaction at 200 keV, These tempe- 
ratxares can be of practical interest in the next generation 
fusion reactors when attairmient of such I'ligh temperatures 
becomes feasible. It is to be noted, however, that for these 
cases, deviation from 3“D Maxwellian leads to a reduction in 
<av>. The same is the case with the results shown in Fig, 3.7 
This has already been explained qualitatively in Section 3.3.1 
in connection with the Tables 3 A to 3.6, 
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Figs, 3.8 to 5,11 show <'ov> vs v.dth 3 as a parameter. 
Since <av> changes fast v/ith T^, a narrov/ range of temp era txxce 
is chosen for each plot, so the difference in <ov> for diffe- 
rent values Ox 3 can be displayed with clarity. The temoera— 
tures chosen are the ones of most practical interest, e,g, 

10 keV for the DT reaction in Fig, 3.8 and ^ 30 keV for 
the DD and DEg^ reactions in Figs. 3.9 and 3.10. It can be 
seen that xor 3=1, the <0v> is lowest, and it is larger for 
both 3 = 0.5 and 2,0. 3 — o yields the highest value of <av> 

followed by the value corresponding to 3 — . Finally, in 
Fig. 3.11, a representative plot of <av> vs is given where 
the reverse is true, 

A suiiLnary of the important conclusions of this chapter 
together with some suggestions is presented briefly in the 
next chapter. 


It 



<<n;> X 10^^ (m^ S'") 



FIG. 3*8 Variation of <cri^>pj Vs with p as a 
parameter. 
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Chapter 4 


SUI''1I4ARY Alffi CUGGESTICNS 


4 . 1 SUMMARY 

After a brief introduction to the present work in 
Chapter 1, we presented a detailed a-nalysis nf the problem 
of calculating fusion reaction rates under non-i'‘Iaxv/e Ilian 
conditions in Chapter 2. 

In Sections 2,1 and 2.2, the basic equations were written 
down and the underlying assumptions and their validity was 
discussed. The basic eq.uations v/ere then used in Section 2,3 
to calcu3-ate <av> for 1-, 2- and 3-dimensional Maxwellian 
distributions. The expressions for these cases are relatively 
easy to establish. To develop a suitable expression for <crv> 
for a general two- tempi eratvire pseudoHyiaxwe Ilian velocity 
distribution (Section 2.5), transformation to CM and relati.-e 
velocity coordinates was done (Section 2,4), The final 
expression was then obtained in Section 2,6, We further 
discussed the basis of comparison for different values of 
<0v> in Section 2,7. In Section 2,8, suitable dimensionless 
parameters were introduced in the various expressions developed. 

In Chapter 3, vai‘ious details of the niimerical computa- 
tions were discussed. The results were presented and discussed 
in Section 3.3. The main conclusions are given below. 
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4.2 CONCLUSIONS 

The calculations reported in this work indicate that 
the fusion reaction rates for the DT, DD and DHe^ reactions 
are enhanced if the velocity distribution deviates from the 
thermodynamic Maxwellian distribution, at temperatures of 
maximum practical interest. For er.ample for the DT reaction, 
the reaction rate may be enlianced by 10 to 20% in the tempe- 
rature range of 5 to 10 keV, Similar conclusion is valid 
for the DKo-^ reaction for 20 to 40 keV range. Enhancement 
in the DD reaction rate is however seen to be rather small 
(see Section 3.3 for details). 

We should stress, however, that the above conclusions 
are based on considering deviations in the form of a two- 
temperature pseudo-Maxwe Ilian distrioution. Their validity 
for arbitrary deviations will need to be verified by addition- 
al computations and/or analytical techniques. 

The calculations reported here also show that at higher 
temperatures (more than 20 keV for DT and more than 40 keV 
for DD and DHe' reactions) , the deviation from Maxwellian 
distr'ibution reduces the reaction rates, in contrast to the 
opposite effect at lower temperature mentioned above. 

4.3 SUGGESTIONS 

1 . In this work we have considered only two-temperature 
pseudo-Maxwe Ilian dis-cxibutions. Arbitrary departures from 
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Maxwsllifln d.is'tr'XloT.i'tj.cn nc'Sd. ■fco tG coiisidGiTGci 'to draw strongGr 
conclusions. 

2, UndGr CGriain condiXioiis (sgg R-si. [“li j) ^ a sirong 
non-Ma;iwe Ilian "tail is superipiposed on "the i-laxwellian distri- 
bution, The effect of this on the fusion reaction rates can 
be investigated very much along the lines of this work. 

3. The validity of the present calculations breaks 
down at the very high densities encountered in inertial 
fusion due to correlation effects. This should be incorpora- 
ted in the conclusions, [see Section 2.2j 

4. As was mentioned in Section 2.2, the polarization 
effects are known to enliance fusion reaction rates by as much 
as a factor of two under suitable conditions. This should 
also be investigated. 

5, For all the calculations v/e have used the data 
shovm in Appendix B, There appears to be some need to improve 
this data in view of the discussions in Section 3.1. Further- 
more, we have used a sharp cut-off at 7 and 1000 keV, This 
range could be enlarged by extending the least squares 
approximation curve. This will have some effect on <av> at 
very low temperatures (1 to 5 keV) or at very high tempera- 
tures (above, say, 100 or 200 koV) . 
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A1 


APPEI’LIy: A 


Proof of dPvy, d^Vo = d''V d-^v [f .e« frojii Eq, 2,4 bj 

i ^ 


From Eqs. 2.16 and 2. 17 5, have, 


V = V — V 
X x^ 


V = V - V 

y Yo Y^ 




Zj z^ 


= '‘a * “b V 


V = d, V -!- V 

y y^ ^ Y'2 


V = d, V, + 
z A z^ £> z^ 


(a1.l) 


where 


'.L 




“A ■'•■ 


and 


nix 


^A ■'■ 


By using xiie definition of Jacooian ; — T 'f’ -i » 'w® 


can wi'ce, 


' ^ (7i ) . 

Iv^ ,V2) 



70 


or 


J = 




j = 





al 




a- 


ai^ 


al., 

^^x 




#X- 

ai 


yi; 

j V 

^ 1 


ai, 

ai, 

^^x 

Zy_ 


11. 

bl„ 





i 

^1 

ai. 

1 



dJx. 



\lt.. 

11. 




i % 





-^12. 



ill 


tv * 

9 yq 


a ^ Vq 

ai2 

^it 

3''x 







■&12 






(al 


The partial derive tiA’-es in 10. c1 ,2 c&.r be substituted 
from Sq. a1.1 to obtain 
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A2 




- d^V d^v 


Proof of -i- = hV^ -:- javf Li.e. Eq. 2.25j 

wher^e i or z 


Vm 


IHr Vo . 
b 2l 


I-n?’ . - 1 - u V 7 
1^ 2. 


(a2.l) 


S_te£s 


From Sqs, 2.16 and 2.17, v;e have 


and 


M 7. = v^. + u,^ 


” ^^21 - ''11 


From v/here we can have \\ . and v,- . as 

1i ii 


^1i = 


IIV. - ITia V. 
J3 _ _ 1. 


and 


I'lV. + HIt:. V. 

■,r _ 1 b 1 

^2i " ‘■‘"'H 


Hence, the left hand side of the a'oove ji-q, a2,1 can be 
written as. 


“A''li ■'■ Vll 


i‘V. - CgV. 2 Vl 2 

% L ‘ yj * ‘ ‘ * b iflp L 


H 
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A. [ m, (h7 - mgv.)^ m3 (KV. m.v.; 

Yr 


\2 




L"'a ■■ “a'^b 


v2 -I- .«aK^v| -;. 


llfvf (m^+mg) + ^“a “b^ •* 


■^r2 


2 

M 7 f ^ V. 


I'l vf + a 


R.K.S. 


Thus, 




"■' Ai 


2 2 
[■: 7-7 + M vr 
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APFEIJDIX B 


VAUJES of and a USED FOR 

DHe^ 

CALCULAIIOEo REPORTED IN THIS ¥ORK (FRON REF.[1 j) 




S.No. ! 
! 

Deuteron Energy, 

E 



i a ,x10^® 

= DHe^ 



keV 

. I 

m2 

1 rr,2 

I { 

I ' 

7.00 

0.00010 I 

0.0 

• 0.0 

I 2. 

f 

8.00 

0.00044 II 

o 

• 

o 

0.0 

3. : i 

f 1 

9.00 

0.00099 i'l 

0.0 

0.0 

4. i ! 

10.00 

0.0018 '! 

0.0 

0.0 

! ' Cj ft 

ii • ' ! 

20.00 

0,055 

0.00054 

0.0 

6, 

30.00 

0.280 i 

0.00222 

0,0 

7. 1 1 

40.00 

0.710 i 

0.005 

0.000245 

8. ' ■ 

50.00 

1.40 i 

0.00849 

0.000924 

9. 

60.00 

2.20 1 

0.013 

0.00225 

10. i 

j 

70,00 

3.05 i 

0,0175 

0.00473 

11. 1 

80.00 

3.85 ‘ 

0.0222 

0.0082 

12. j 

90.00 

4.45 I 

0.0270 

0.0135 

= 13. • 

100.00 

4.80 ■ 

J 

0.0315 

0.0197 

■14. 

200.00 

4.00 j 

0,0640 

0,185 

15. 1 

1 

300.00 

2.30 

0.0848 

0.480 

He, 

400.00 

1.40 

0.108 

0.724 

■17. j 

« 

500.00 

1.05 i 

I 

0,120 

0.780 

:i8. 

600.00 

0.750 1 

0.135 

0.648 

19. 

700.00 

1 

0,580 

0.141 

0.490 

20. 

j 800.00 

1 0.450 

0.152 

0.369 

21. 

1 900.00 

1 0.350 

0.170 i 

0.300 

,22. 

1000.00 

i 0.510 

\ 

J 

0.180 ! 
1 1 

0.244 
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At’i’EMLIX 


C 


TABLE OF SOLE COl^oTANTS AND 


CONVERSION FACTORS 


USED IN THIS WOR[^ 


nij^ -- 2.0141 u 

= 3.01605 u 

111:, = 3.0160^ u 
Ile^ 

1 u = 1.660438 X 10"^'^ Kg. 
1 keO = 1.602 X 10“'^^ u 


h = 

6, 

.626 

X 

10 j-a 

N = 

0 

6. 

.022 

X 

10^^ 

e = 

1. 

.602 

X 

10“"'^ C 

^0 = 

8, 

.854 

X 

a 

r- 

OJ 

0 

CM 

1 

0 

C = 

2. 

,998 

X 

10® ms ^ 


1 barn = 10 




7i' 


[1J 


[ 2 ] 


[3J 

[4] 



[6J 

[7] 

[ 8 ] 
[9J 
[ 10 ] 
[ 11 ] 


[ 12 ] 

[ 13 ] 
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